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ABSTRACT 

We present the first X-ray analysis of the diffuse hot ionized gas and the point 
sources in IC131, after NGC604 the second most X-ray luminous giant H II region in 
M33. The X-ray emission is detected only in the south eastern part of IC131 (named 
IC131-se) and is limited to an elliptical region of ~200 pc in extent. This region ap- 
pears to be confined towards the west by a hemispherical shell of warm ionized gas 
and only fills about half that volume. Although the corresponding X-ray spectrum 
has 1215 counts, it cannot conclusively be told whether the extended X-ray emis- 
sion is thermal, non-thermal, or a combination of both. A thermal plasma model of 
kT e = 4.3 keV or a single power law of T ~ 2.1 fit the spectrum equally well. If the 
spectrum is purely thermal (non-thermal), the total unabsorbed X-ray luminosity in 
the 0.35 - 8 keV energy band amounts to Lx = 6.8 (8.7) x 10 35 erg s _1 . Among other 
known Hn regions IC131-se seems to be extreme regarding the combination of its 
large extent of the X-ray plasma, the lack of massive O stars, its unusually high elec- 
tron temperature (if thermal), and the large fraction of Lx emitted above 2 keV (~40 - 
53%). A thermal plasma of ~4keV poses serious challenges to theoretical models, 
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as it is not clear how high electron temperatures can be produced in H II regions in 
view of mass-proportional and collisionless heating. If the gas is non-thermal or has 
non-thermal contributions, synchrotron emission would clearly dominate over inverse 
Compton emission. It is not clear if the same mechanisms which create non-thermal 
X-rays or accelerate CRs in SNRs can be applied to much larger scales of 200 pc. In 
both cases the existing theoretical models for giant H II regions and superbubbles do 
not explain the hardness and extent of the X-ray emission in IC131-se. We also detect 
a variable source candidate in IC131. It seems that this object is a high mass X-ray 
binary whose optical counterpart is a B2-type star with a mass of ~9M Q . 

Subject headings: ISM: H n-regions — ISM: superbubbles — galaxies: individual 
(M33) — X-rays: individual (IC131) — X-rays: ISM 



Introduction 



This study is part of the Chandra ACIS Survey of M33 (ChASeM33, IPlucinskv et all 120081) 
and presents the first deep, high resolution X-ray images of IC 1 3 1 . Within ChASeM33's sensitivit y 
limit of ~2 x 10 34 ergs" 1 (0.35- 8 keV) IC131 is, along with NGC604 dTullmann et al. Il2008h . 
the only giant H II region (GHR) in M33 which shows significant diffuse extended X-ray emission. 



Because of their large sizes, masses, stellar contents, and luminosities, GHRs (cf. iKennicutt 



19841) provide excellent laboratories for studying the relation between stellar feed back and the 



evolution of the indiv idual com ponents of the IS M inside them. GHRs like 30 Dor (TWang II 1999c 



Townslev et al. 1120060 . N51D dCooper et al.ll2004f) . or NGC604 dTullmann et al. 112008b are strong 



X-ray emitters thanks to their massive O star population. The heavy mass loss of these stars can 
produce strong colliding winds which, together with contributions from SNe and SNRs, shock- 
heat the gas and forces it to emit thermal X-rays at temperatures around kT e = 0.6 ± 0.2keV. In 



to a thermal component, also a non-thermal component (e.g. 


. Wolk et al. 


|2002 


^ — 1 — , --- — 
: Bamba et al. 


2004; 


Cooper et al. 


2004; 


Smith & Wang 


2004; 


Muno et alJ2006: 1 


vladdox et alJ2009b. The non-thermal 



radiation could be explained b y synchrotron emission, inverse Compton scattering, or by particle 
acceleration in shock regions dParizot et al J 12004c iButt & Bykovll2008l) . In all cases, the X-ray 
emission could be linked to massive stars or their successors. The issue one faces with GHRs is to 
find a mechanism which can ionize the gas on spatial scales which can easily exceed that of normal 
H II regions by one to two orders of magnitude and is consistent with the stellar energy feedback, 
X-ray luminosity, gas mass, and other predictable quantities. 



Although IC131 has been observed in numerous M33 surveys across almost all wavelengths 
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Rosa et all 1984 


: Landsmanetal.il 1992 


:Lons et al. II 1996: 


HiDDelein et alJ2003 


Pietsch et al. 


2004 


: iMassev et al. . 


2006; 


Tabatabaei et al. 


2007; 


Plucinskv et al. 


2008 


), little is known about the 



stellar population, its age, and the different components of the ISM. The ChASeM33 data set, sup- 
plemented by other observations at optical and infrared wavelengths, is the first data set that has 
the combination of S/N and spatial resolution to permit a detailed study of the X-ray emitting gas 
in IC131. Therefore, the primary focus of this work is to study the hot ionized medium (HIM) by 
constraining basic plasma parameters, to determine the likely ionization mechanism of the hot gas, 
and to compare properties of the gas to other star forming regions to learn more about the evolution 
of the gas in GHRs. The second purpose of this study focuses on the analysis of the point source 
population in IC131. 

In the following, we give an overview of the observations used in this study and the data 
reduction steps that have been applied (Section 0). We then present high resolution (2") images of 
the diffuse X-ray emission in different energy bands and compare these to multi-wavelength data 
to constrain the morphology of IC131 (Section [3TTb . A spectral analysis is carried out in Section 
I3.2l while Section [3~3l discusses the possible origin and ionization mechanism of the HIM. Section 
@]deals with the analysis of the X-ray emitting point sources detected in IC131 and Section[5]gives 
a concluding summary of the most important results. 



2. Observations and Data Reduction 



We have chosen for our analysis 3 different ACIS-I fields, namely Fields 3 and 4 from ChASeM33 
(consisting of ObsIDs 6380, 6381, 6382 6383, and 7226) and one archival field (ObsID 1730). 
The location of each field is displayed in lPlucinsky et al.1 (|2008|) . In the ChASeM33 observations 
IC131 is approximately 4'.3 off-axis and in ObsID 1730 ~9'. All observations were performed 
with ACIS-I in 'Very Faint' (VFAINT) mode and all data were reprocessed with CIAO 4.0.1 and 
CALDB version 3.4.5 . Basic data reduction was carried out according to the steps outlined in 
Tullmann et al. I d2008|) . The GTI-corrected total exposure time is 422 ks. 



We constructed exposure-corrected X-ray images from the merged event list in the follow- 
ing energy bands: 0.35-8.0keV (broad), 0.35- 1.1 keV (soft), 1.1 -2.6keV (medium), and 2.6- 
8.0 keV (hard). All images were binned by 4 pixels to a resolution of ~2" (or 7 pc) to reduce noise 
fluctuations in the binned pixels. Since IC131 was observed at off-axis angles >4', the radius of 
the PSF (90% encircled energy at 1.49keV) in our images is between 2 / .'l and 8". Hence, binning 
to a pixel size of 2" x 2"does not degrade the angular resolution significantly. 



In addition to the X-ray data, we also employed archival broadband (R , V) and narrowban d 
imagery (Ha, [O in] A5007A) from the Local Group Galaxies Survey (LGGS. lMassev et al. II2006I) . 
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the MMT (r, g, u), and Spitzer (IRAC/MIPS images at 3.6, 8, and 24/mi) which allow a multi- 
wavelength comparison of the different phases of the ISM. Because we use these images merely 
for a morphological comparison, only standard data calibration was applied, except for the Ha and 
[O ill] images, which have R and V band continuum images subtracted. 

We extracted spectra from the diffuse X-ray emission after point source removal. Because 
the IC131 pointings do not share common roll and pointing directions we did not create spectra 
from the merged event list. Instead we extracted source and background spectra for each i ndivid - 
ual event list and merged them by following again the steps described in ITullmann et al. I (|2QQ8f> . 
The resulting source spe ctra were backgro und subtracted and then grouped using the grouping 
scheme of ACIS Extract (IBroos et al.ll2002l) until a signal to noi se ratio (S/N) of 3 per channel was 
achieved and fitted with XSPEC ver . 12.4.0 (Arnaud Ill996|) . Uncertainties are given on a 90% 
confidence level and were calculated with the error routine. 



The spectrum of the only X-ray detected point source in IC131-se, CXO J013315. 10+304453.0 
(here after also FL073 as it is source number 73 in the ChASeM33 First Look catalog of lPlucinsky et al. 
20081) . was extracted in a similar manner. This source is located at RA 01 h 33 m 15!10 and Dec 
30 o 44'53'.'00 (see red cross in Fig.[T^) at the western tip of the extended X-ray emission. 
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Fig. 1. — (a): Continuum-subtracted Ha image of IC131 from the LGGS. (b) Same as a), but for 
[O m]A5007A. (c): Spitzer IRAC and MIPS data obtained at 3.6/xm (red), 8/xm (green), and 24/zm (blue). 
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(d): Multicolor image, showing Ha in red, [0 m]A5007A in green, and X-rays (0.35 -8.0 keV) in blue, (e): 
High resolution (2") Chandra X-ray images (red: 0.35- 1.1 keV, green: 1.1 -2.6 keV, blue: 2.6-8.0 keV). 
(f): Optical false color composite consisting of MMT broadband r, g, and u imagery. Positions of three 
spectroscopically classified WR stars (see labels in panel f) for their spectral type) and 21 unclassified stars 
within the bubble (yellow circles) are taken from Ma ssey et al. zoom into the region surrounded 

by the yellow dashed box is shown in Fig. [5] (see Sect. 4.1) and reveals a possible optical counterpart of 
FL073, the bright X-ray point source marked by a red cross in panel e). Contours in panels d) -f) represent 
IR emission at 3.6^m (red, 0.42 - 1 .50 mJy/sr) and 24/xm (white, 3 1 .9 - 39 mJy/sr). See the electronic edition 
of the Journal for a color version of this figure. 



3. Diffuse emission 

3.1. Imaging Results 

From the Ha and [O m]A5007A imaging shown in Figs.[T^, b, and, d one can see that IC131 
consists of two prominent regions, a hemisphere in the north west and a hemisphere in the south 
east. The brightest optica l emission in the form er region seems to be related to an H II region which 



is ionized by a Bl star (Ma ssey et al. 1 120061) . whereas the latter complex is the most dominant 
feature in the optical. This region appears to be a superbubble of ionized gas with a radius of 
~100pc that opens to the SE (called IC131-se, hereafter). These two complexes appear to be 
separate bubbles and give IC131 a characteristic dumbbell or hourglass shape. The western edge 
of IC131-se shows enhanced Ha and [O ill] emission, running roughly from RA 01 h 33 m 15!0 and 
Dec 30 o 45'25'.'0 to RA 01 h 33 m 14 f 50 and Dec 30°44'30"0. This area can be associated with young 
stars, among them are three WR stars, which seem to have formed along the edge (see yellow 
diamonds in Fig. QJ for their positions and type). It appears likely that these stars photoionize the 
gas and cause the optical enhancement. Moreover, the IR emission at 24 fim (white contours in 
Figs, [lji -f), a tracer of warm dust, is well aligned with the western edge of the bubble and seems 
to follow the stellar distribution. From Figs. [Hi and e it can be seen that except in the east, the 
diffuse X-ray emission seems to be confined by the superbubble and that it extends eastward, away 
from the loci of the hot stars. A peculiarity in the large scale distribution of the X-ray emission 
in IC131-se is noteworthy. From the image of the 0.35-8.0keV emission shown in Fig. [T]i, it 
appears as if the diffuse gas is confined to the southern half of the superbubble. Such a distribution 
is unexpected, especially if the gas was produced by stellar winds, SNRs, and SNe, because the 
high-pressure and supersonic ally streaming gas should fill the bubble almost completely. Shearing 



by differential rotation (e.g., iPalous et alj|1990|) is an unlikely explanation, as that would equally 
affect the bubble and the HIM and is likely to act perpendicular to the major axis of the HIM. It 
would be very illuminating to know how such a lopsided morphology can be produced. 
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Using funcnts from the FITS utility package FuntoolsQ, which counts the number of 
photons in a specified regions, we determined that ~85% of the detected diffuse X-ray photons 
in IC131-se have energies >lkeV (see the greenish emission in Fig. and 34% are emitted 
above 2 keV. For NGC604, the second largest H II region in the Local Group, the corresponding 
values are 45% and 10%, respectively, and clearly indicate that the X-ray emission in IC131-se is 
significantly harder. 

Although there are some NIR point sources in IC131-se seen superposed on the X-ray emis- 
sion, the extended NIR emission and the X-rays are anticorrelated, suggesting that they do not 
occupy the same volume. This anticorrelation cannot easily be explained in terms of absorption, 
since most of the X-rays we see have energies in excess of 1 keV and should not be significantly 
absorbed. 

Among the 21 identified stars located in side the bubble (yel low circles in Fig.Q]), none were 



spectroscopically classified as type WR or O (Ma ssey et al. IfeOOq) . This implies that most massive 



stars inside the bubble have already exploded and that we face IC131-se in a relatively late evolu- 
tionary stage, i.e. close to the end of the life cycle of massive stars. The SNe events could have 
triggered the formation of a sparse, new generation of stars which see med to have formed along 



the western edge of the bubble; among them are three WN-type stars (IMassey et al. Il2006|) . This 
new star formation site is most likely associated with an extended reservoir of molecular gas which 
could be pristine gas from the remainder of the parent molecular cloud. It appears as if star for- 
mation progresses from east to west and has already consumed a large part of the molecular cloud 
(provided the bubble formed from the same m aterial). In this picture, SNe shocks or stellar winds 



as well as the hot gas itself (|Smith et al.111996 ) could have swept away or destroyed any significant 
preexisting accumulation of dust from the bubble's interior, which would also explain the lack of 
IR emission. 

Apparently, IC131-se is very different from 30 Dor, NGC604 or any other known superbubble 
in at least two important aspects: first, IC131-se does not show a clear spatial correlation between 
young stars and diffuse X-ray emission and second, the X-ray emission in IC131-se (regardless 
whether thermal or non-thermal) seems to be much harder on spatial scales that have not been 
observed in other H II regions or superbubbles before. The first aspect is explainable by different 
ages of the stellar population, the latter, however, remains puzzling and requires an analysis of 
possible emission mechanisms. 



'httpV/www.cfa.harvard.edu/^john/funtools/ 
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Fig. 2. — In the left panel the X-ray spectrum of the diffuse emission is fitted with a soft (dotted 
line) and hard (dashed line) thermal plasma model, whereas the right panel shows the best-fit of a 
non-thermal power law model. 



3.2. Spectroscopic Results 

As discussed in the previous section, the X-ray spectrum of the superbubble in IC131 is much 
harder than usually observed in other typical GHRs, such as NGC604. To explore this, we extracted 
a spectrum of the diffuse emission, after excluding a region around the point source, and performed 
spectral fits to help us understand the nature of the X-ray emission. This spectrum contains 1215 
counts and is shown in Fig. [Altogether with the best-fit models for a purely thermal plasma (albeit 
two thermal components) and a non-thermal plasma (see also Table [I])- Although the spectrum is 
relatively smooth, there are some features at soft energies around 0.68 and 0.9 keV which suggest 
the presence of O vill, Ne IX, or Fe L emission. At energies above 2 keV the spectrum looks 
relatively smooth and there is no real indication of emission lines, such as Fe Ka at 6.7 keV. 
However, it should be kept in mind that Chandra's sensitivity is low at these energies. 

By adopting an alternative grouping scheme provided by CI AO (dmgroup, NUM_CTS=20), 
we tried to increase the number of groups in order to improve the spectral resolution around the 
energies of the potential emission features (see Fig. [3]). The resultin g spectrum was fi tted with 



a modified apec model which has continuum, but no line emission (|Smith et al.ll2001l) and two 
unconstrained Gaussian profiles at 0.68 and 0.89keV. Because this model fits the emission line 
features well and because the line widths of the Gaussians are small (< 0.02 keV), the presence of 
a thermal component seems to be indicated. However, we also fitted the spectrum shown in Fig. [3] 
with an absorbed power law and found that the residuals are barely significant (Fig. 4). Although 
the present data are suggestive of a thermal component, it is not required. 
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Table 1 . Fit parameters (diffuse emission) 



tbabs x (tbvarabs x . . . ) a 
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4.39^| 7 
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h0.09 
-0.21 
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a The tbabs model 


assumes a 


fixed 


Galactic Nh column density of 6 x 10 20 cm 2 



(jDickey & Lockman Ill990h and solar abundances, whereas the tbvarabs model accounts for 
internal absorption in M33 and uses abundances of Z = 0.6Z Q . b Under the assumption that the 
model is a good fit to the data, the P value is the probability of getting a x 2 which is greater than 
the observed x 2 bs- c Ki 2 are the normalization constants of the fit for the individual model com- 
ponents. ^Photon index. e kT 12 represent the electron temperatures of the individual model com- 
ponents. ^Absorbed (F x ) and unabsorbed (F x ) fluxes are tabulated in the 0.35- 8.0 keV energy 
band. ^Unabsorbed luminosities in the 0.3 5 - 8.0 keV and 2.0 - 8 .0 keV energy band, respectively, 
assume a distance to M33 of D = 817kpc (|Freedman et al. 11200 1|) . 



- 10- 



grouped NUM_CTS=20, tbabs(tbvarabs(apec_noline+gau+gau)) 

















++++++ i+ttt^i 



0.5 



Energy (keV) 



Fig. 3. — This plot shows part of the diffuse spectrum which was grouped with dmgroup until 
each group had at least 20 counts. It was fitted with an absorbed thermal apec model from which 
the line emission was removed (dotted line) and two Gaussian profiles (dashed lines). 



grouped NUM_CTS=20, tbabs(tbvarabs(pow)) 
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Fig. 4. — If we fit a simple absorbed power law to the same spectrum shown in Fig.[3l the residuals 
at 0.68 and 0.89 keV do not appear significant enough to claim the existence of emission lines, 
implying that the power law is an acceptable fit, too. 



Since the line features (if they exist) are weak, we also fitted the spectrum with a number of 
simple models using thermal and non-thermal components. The results of the best model fits are 
listed in Table [D In all cases the Galactic H I col umn density was modeled with the tbabs model 
(|Wilms et al.ll2000|) and fixed at 6 x 10 20 cm -2 (|Dickey & Lockman Ill990|) . For the absorption i n 
M33, we used the tbvarabs model, adopted a metal abundance of 0.6Z & (IVilchez et al. 1988). 
and set the N H initially to 1 x 10 21 cm~ 2 to account for the average value in M33 (INewtonlll980|) . 
The iVn in M33 was allowed to vary freely during the fitting process. For the apec model, the 
metal abundance was also set to 0.6Z Q . 

From a statistical point of view all of the models shown in Table[T]are acceptable. An absorbed 
thermal plasma model appears to be as likely as an absorbed power law. The predicted electron 
temperature of the thermal model of about 4.3 keV would be the highest determined so far for any 
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GHRs. However, the unabsorbed luminosity of ~7.0 x 10 35 erg s -1 is about half that of NGC604. 
About 50% of the absorption-corrected X-ray luminosity of IC131-se is emitted in the 2-8keV 
energy band, which indicates that the emission is much harder than in "nor mal" giant H II regions , 



such as NGC604 or 3 Dor, for which the corresponding fraction is <5% (|Townsley et al. 



2006; 



Tullmann et al. 1 120081) . As there is no physical reason to believe that the thermal plasma within 
the complex environment of a GHR would have a single temperature (1-T), we also considered 
a two-component thermal model (2-T). The best fit has a soft (0.14keV) and a hard (3.4 keV) 
component and is, not surprisingly, somewhat better than the one from the single temperature 
model. This second component allows for a significantly better fit of the soft emission below 
0.8 keV and likely yields a more accurate estimate of the X-ray luminosity. The temperature of the 
higher temperature component in the 2-T model is lower than the temperature for the 1-T model, 
which is likely a consequence of the larger N H of the 2-T model. However, the temperature of the 
higher temperature component of the 2-T model is consistent (within errors) with the temperature 
of the 1-T temperature model. We also fitted the spectrum with NEI and bremsstrahlung models, 
which also produced acceptable fits. However, they are not well motivated physically, given that the 
ionization timescale of about 4 x 10 13 s cm~ 3 as derived from the former model implies ionization 
equilibrium conditions and the latter considers the unrealistic case of a collisional plasma that is 
pure H and He. 

Non-thermal X-rays, either through synchrotron emission or inverse Compton scattering are 
well known from SNRs and H II regions. This, together with the fact that the spectrum of IC131-se 
appears relatively smooth above 2 keV, motivated us to also fit the spectrum with a purely non- 
thermal model and a combination of thermal and non-thermal plasma models. Like a single apec 
model, the single power law provides a good fit to the data, too. The photon index of T = 2.13^q ^ 
is within a reasonable range expected for SNRs and superbubbles with non-thermal X-ray emission 
and the is consistent with the integrated line-of-sight column density. 

We also considered combinations of thermal and non-thermal emission, which seemed in 
some sense natural, given that a combination of thermal and non-thermal emission has also been 
reported for other GHRs and superbubbles. These models also led to statistically acceptable 
fits, and the best fit has a photon index of T = 2.66 and a very soft plasma temperature of 
0.07 keV. However, the best-fit column density of iV H = 8.4 x 10 21 cm~ 2 appears unphysical. 
Bal mer line extinction measu rements in IC131-se indicate that reddening is low and mostly Galac- 
tic (Relano & Kennicutt 1 l2009h . Likely as a consequence of the large column density in this model, 



the unabsorbed luminosity in the broad energy band is unacceptably high, about a factor 60 higher 
than predicted by the other models and almost a factor of 30 larger than in NGC604. Constraining 
iV H to acceptable values tends to drive the best fit to a simple power law with little or no contribu- 
tion from the thermal component, which is understandable since the pure power law fits actually 
do fairly well in fitting the overall shape of the spectrum. We also tried a combination of an apec 
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and a srcut ([Reynolds & Keohandll999|) model. This best-fit model predicts a luminosity which 
is comparable to that of 30 Dor and a spectral radio index which is much higher than the observed 
one of T = 0.3 (F. Tabatabaei, priv. comm.). Given the limited number of counts in our spectrum 
and the fact that all of the models we considered are very simple, we do not believe that the combi- 
nation of thermal and non-thermal emission should be ignored, although the model results do not 
support that picture. 

To ensure that the emission was not contaminated by background sources or insignificant X- 
ray point-like sources, we also removed the X-ray emission within the cone which coincided with 
the IRAC 3.6/im emission and refitted the spectrum. No significant differences to the previous best 
fit results were obtained. Obviously, one would like to obtain a better X-ray spectrum of IC131- 
se, but that may be unrealistic for the foreseeable future, because our Chandra exposure totaled 
422 ks. Unfortunately, the available XMM-Newton observations provide no further insights as the 
exposure time is either too low to extract a spectrum, IC13 1-se falls onto a chip gap or is more than 
8' off-axis, where it is strongly vignetted. 



3.3. Discussion 



For a GHR or superbubble the X-ray properties of IC 13 1-se are extreme. The fraction of 
the hard X-ray energy flux above 2keV exceeds that of 30 Dor and NGC604 by one order of 
magnitude (50% vs. 5%). If the emission is thermal, the temperature of the hot gas in IC 131- 
se is 4.3 keV, at least a factor o f seven higher comp ared to t hese and other kno wn Hn regions 
and superbubbles, such as Nil dMaddox et alJlioOa) . N51D Jcooper et alJb004h . and 30 Dor C 
(Ba mba et al. I l2004l : lYamaguchi et al.ll2009l) . This raises the question whether the X-ray emission is 
a relic of the same process that caused the expansion of the bubble or if it is due to processes which 
happened in the recent past. In order to answer that and to constrain the nature of the emission, 
we evaluate the possibilities that the emission is either thermal, non-thermal, or a combination of 
both. Two other factors need to be incorporated into any answer to this question. IC131-se and the 
X-ray emitting region are large, a little more than 200 pc across, and this region contains few if any 
O stars to provide a continuing source of energy to power the emission. 



3.3.1. Is the X-ray emission thermal? 

There are many reasons to believe that the X-ray emission from IC131-se is due to emission 
from a thermal plasma. Firstly, the X-ray emission from most GHRs is dominated by emission 
from thermal plasmas. Secondly, a well-developed explanation for the X-ray emission exists in 
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the form of shocks created by a combination of stellar winds and SNe explosions. Thirdly, the 
X-ray luminosity of IC131-se is not extreme and is comparable to other superbubbles with thermal 
emission. Fourthly, from an observational perspective, the single absorbed apec model with a kT e 
of 4.3 keV provides an adequate fit to the data, and of all of the thermal models which do so, a low 
value of ]Vh is predicted that is consistent with other observations. Furthermore, there are hints of 
line features at positions expected from a thermal plasma. 

The accepted interpretation for X-ray emission from Hll regions and superbubbles (e.g. , 



Weaver et al. 1 1 19771 : McCrav & Kafatos 1119871: MacLow & McCrav 1988; Shul l & Sakenlll995h 



is that it arises from a thermal plasma that has been heated by supersonic colliding winds from 
massive stars, SNe, and SNRs to temperatures <1 keV. Within this context the data discussed in 
Section 1341 suggest that in the case of IC131 a large superbubble formed by the combined action 
from stellar winds, SNRs, and SNe, filling its interior with a hot X-ray emitting gas, similar to 
30 Dor or NGC 604. If this interpretation is correct, the absence of a hot and massive stellar pop- 
ulation within the bubble implies th at these stars already exp loded as SNe and the superbubble is 



a relatively old structure (>8 Myr. LMaeder & Meynet II 19 94). It might also suggest that the south 



eastern part of the superbubble was swept away by SNe blast waves. Alternatively, the bubble 
could have become unstable after the expansion slowed down and cooling set in, which finally led 
to large-scale breakout of hot gas into the ambient ISM. Either case suggests that IC131-se is in an 
advanced evolutionary stage. This is different from NGC604, the other X-ray bright GHR in M33, 
because in this case the present stellar generation in the western part appears to be too young to 
produce a significant number of SNe. An evolved and aged superbubble would also explain why 
no SNRs have been reported for IC131-se, they could have been thermalized a long time ago. It 
could also be that recent SNe have escaped detection, because they are evolving into a low density 
medium. 

The main difficulty with a thermal interpretation of the emission from IC131 is the high 
electron temperature. Hot X-ray plasmas (kT e > 3 keV) extending on scales of ~200pc have not 
been reported for GHRs so far. If the plasma is purely thermal, IC131 is the first such region, 
allowing us to study gas under extreme conditions and to learn more about the mechanism, which 
can produce such hot plasmas. Before attempting to answer whether the conditions under which 
the plasma might be heated to these temperatures exist, it is worthwhile to explore whether such a 
hot plasma is allowed based on energy or mass arguments. 

The electron density can be calculated from the normalization constant of the best-fit model 
(here the single temperature apec model) and amounts to n e = (0.045 ± 0.014) f^^ 2 cm" 3 (with 
fx being the filling factor). This value is very low (a factor ~5 lower compared to NGC604) and 
would also support the idea that blowout in IC131-se occurred. 
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The internal thermal energy within the superbubble can be written as: 



E th = 2.88n e kT e Vl 



X; 



(1) 



where n e is the electron density, k is the Boltzmann constant, T e is the assumed electron tempera- 
ture of T e = (5.0 ± 1.3) x 10 7 K, and V is the volume occupied by the HIM, which we adopt to be 
(7.0 ± 1.0) x 10 61 cm 3 (an ellipsoid with half-axis radii of 1 lOpc x 72 pc x 72 pc with 50% uncer- 
tainty in each axis). The internal energy is calculated for a fully ionized plasma of Hydrogen with 
10% Helium and a total number particle density of n ~ 1.92n e . We estimate the internal energy of 
the bubble to be (5.9 ± 2.9) x 10 52 erg, assuming a filling factor of fx = 0.9. Because the average 
conversion efficiency between the kinetic and thermal energy in a SN explo sion is unlikely to b e 
one, but closer to 0.5 during the first few 10 3 yr of the evolution of a SNR (ITang & Wan g 2005), 
the thermal energy content in IC 1 3 1 - se should roughly correspond to 118±58SNe, as suming each 
SN releases a kinetic energy of 10 51 erg. Is the number of SNe also consistent with the number of 
O-stars expected from the initial mass function (IMF)? 

A stellar population analysis of the stars along the line of sight to IC13 1-se does not currently 
exist, but the available optical photometry of the stellar population inside the bubble can provide 
some clues on the nature of the X-ray emission. Among the 21 stars seen in projection against the 
bubble, five appear to be located at the edge of the new star formation site which coincides with the 
eastern edge of the molecular cloud. Because these stars are possibly not members of the stellar 
population ins i de the bubble, we excluded them from the following calculations. According to 
Massey et al. I (|2006|) . the remaining 16 stars have apparent (absolute) visual magnitudes ranging 
from 20. 1 3 (-4.65) to 22.99 (-1.78) which implies that their spectral types range from B 1 III to B6 III 
(7< M/M^ < 20) (ISchmidt-Kalen Il982|) and all stars, except one, are potential SN progenitors 
<rWooslevetalJl2002 ). 



These findings confirm our previous consideration that we observe IC 13 1-se during an ad- 
vanced evolutionary stage in which all O stars already exploded and SNe and SNRs heated the 
gas. Based on the available optical data we can narrow down the age of the stellar population 
inside the superbubble. Making use of the stellar mass-luminosity relation and assuming that the 
least massive B star has a mass of ~7 M and that the lifetime of the least massive O star is about 



8Myr (|Maeder & Meynet 1 11994). the age of the stellar population ranges fro m at least 8 M yr to 
77 Myr. By adopting a Salpeter IMF with slopes ranging from 1.6 < a < 3.0 (IKroupa!l200l|) . we 
determine the number of expected O stars with masses 10 < M/M Q < 120 to range from 9 to 21 
(it needs to be kept in mind that the normalization constant of the IMF may seriously affect the 
expected number of O stars if the observed sample suffers from incompleteness). This estimate 
is inconsistent with the number of SNe derived from the internal energy content of IC131-se and 
argues against a t hermal origin of the X-ray emission. Compared to NGC604, which hosts about 
200 O-type stars (|Hunter et al.l ll996) and is ~1.5 times larger, the number of SNe expected for 
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IC131-se is ~130, which is considerably larger than the estimate from the IMF. So either the IMF 
estimate suffers from incompleteness and the number of O stars was significantly larger or the 
diffuse X-ray emission is not predominantly thermal. 

The energy density of the gas amounts to ~8.0 x lCT 10 erg cm" 3 and the corresponding 
cooling timescale is r coo i ~ 3.5 kT e /A(T) n e , where A(T) is the normali zed cooling function, 
whic h amounts to A(T) ~ 1.7 x 10~ 23 erg cm 3 s" 1 for half solar metallicities ([Sutherland & Dopita 
1993b . Hence, the cooling timescale is ~1 Gyr and much greater than the age of 8 -77 Myr of the 
superbubble. This implies that, due to the low density, the plasma remains hot for a long time and 
cannot cool efficiently. However, the main problem remains to be solved, namely how to get the 
plasma heated up to electron temperatures of about 4 keV. 

The total gas mass that one calculates based on a thermal model is reasonable. The ionized 
gas mass traced by X-rays can be written as: 



M x = 1.17 m H n e V7x, 



(2) 



with 1.17 x m H being the hydrogen mass per electron, accounting for the contribution of 10% He 
and assuming fully ionized gas of solar abundances. Adopting again a filling factor of the HIM 
of 90%, the total X-ray gas mass amounts to M x ~ 3300 ± 550M o , which is about 50% of that 
in NGC604. This is considerably greater than the mass of the progenitors of the 118±59 O-stars 
needed to heat the plasma, but conductive mass loss from the walls of the supershell can easily 
account for the mass deficit. 

Hot plas mas are expec ted to expand rapidly and eventually break out from their host envi- 
ronment (e.g.. lTanakall2002|) . The expansion of the X-ray gas is clearly suggested by a compar- 
ison between the internal pressure of the bubble and the external pressure from the surrounding 
ISM. The thermal pressure of the hot gas inside the bubble can be estimated to be P t h/k = 
1.92 n e T e = 4.5 x 10 6 Kcm~ 3 , assuming a filling factor of 90%. This pressure is more than a 
factor of 10 -20 higher than in the Nil and N51D superbubbles and almost 2 times higher than 
in NGC604 (ITullmann et al. Il2008|) . Compared to the canonical value for the external pressure of 
P e / k = 10 4 K cm" 3 of the ISM, a rapid expansion of the HIM is expected and even breakout seems 
possible given the non-detection of cooler gas which confines the HIM towards the east of the bub- 
ble. A bubble which suffers breakout would also be consistent with an evolved age of IC131-se. If 
we assume that the HIM broke out from the bubble, it remains to be understood why the internal 
pressure is still so high. It could be that we fortuitously witness breakout at a very early stage where 
the pressure is still high. Assuming that this phase is short compared to the bubble's lifetime, it 
appears rather unlikely that we observe such a rare event. In case we do, it could explain why we 
have not observed similar objects before. Although a thermal origin is expected from the standard 
H II region or superbubbl e models, electron t emperatures of 4keV are clearly not. Even relatively 
young SNRs like E0102 faughes et allboooh or 0509-67.5 Jchen et alJIioOS) have spectra which, 
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when characterized by simple thermal models, have apparent electron temperatures significantly 
below 1 keV. It is in particular hard to understand how one could create a very hot plasma in a 
situation where a substantial portion of the gas mass was conducted from the walls of the region. 
Therefore, the fundamental theoretical challenge with a thermal explanation is to understand how 
electron temperatures of the observed order can be obtained. Shock waves in SNRs can in principle 
produce ion plasma temperatures between 10 7 - 10 8 K (860 -8600 eV) in regions behind the reverse 
shock. However, the electron temperature should be about three orders of magnitude lower due to 
mass-proportional heating, which is too low to produce X-rays. So any SNR seen in X-rays must 
have some amount of collisionless heating. Two such examples, for w hich the heating app ears 
not to be p roportional to the ion to electron mass ratio are SN 1006 (e.g.. [Laming et al.lll996l) and 
DEML71 ( Rakowski et al. 2009 ). In the se cases the Buneman instability (|Bunemanl 1958 ) and the 
excitation of hybrid waves jLesch 1990 ) in non-radiative shocks are assumed to raise the electron 
temperature by about a factor of 300. It remains unclear whether these processes also occur in 
IC131-se and questionable if the additional heating is large enough to produce the observed high 
electron temperatures. 

If the age of IC131-se is long compared to the Coulomb equilibration timescale, the electron 
temperature should be high, as the protons would have had enough time to transfer their energy 
to the electron s via Coulomb collisions. The estimated equilibration time t eq = 6.75 x T^ 2 /n e 



(ISpitzerl 1 1 97 8[) for IC131-se amounts to 1.5 Myr. This time is much shorter than the estimated 
age of 8 -77 Myr and suggests that the high temperatures are caused by Coulomb heating of the 
electrons. However, the Coulomb timescale in NGC604 is ~ 0.025Myr and sub stantially shorter 
than the age of the stellar population of ~3 Myr ( Gonzalez Delgado & PerezlboOO ). yet the electron 
temperature there is much lower than the one predicted for IC131-se by the thermal model. Thus, 
it remains puzzling why the electron temperature would be so much higher in IC131-se than in 
other GHRs. 



3.3.2. Is the X-ray emission non-thermal? 



Based on the fits, a pure non-thermal plasma is statistically as justified as a pure thermal 
plasma to explain the origin of the X-ray emission in IC131-se. Non-thermal X-ray emission is 
observed in increasing numbers of SNRs, ranging from youn g objects like Cas A and SN1006 (e.g., 
Koyama et al. 1995 : Long et al. 20031: iHelder & Vink 2008 ) to older objects with slower shocks, 
such as G156.2+5.7 iKatsuda et al.ll2009h. In such objects, particles appear to be accelerated by 
diffusive shock acceleration ( Bell ll97"8TMalkov & O'C Drury 2001 ) and non-thermal X-rays are 
most likely due to synchrotron emission or inverse Compton scattering. 



Because SN1006 or 30 Dor C show distinct cap-like regions which emit strong synchrotron 
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X-rays, one could therefore speculate that the hard X-ray emission in IC13 1-se is similar in nature, 
except that in IC131-se the emission is projected towards the observer so that it is seen end-on 
rather than edge-on. Non-thermal emission has also been invoked t o explain the hard tails in the 
spect ra of other H II regions and super bubbles in the Milky Way (e.g jMuno et all2006 : Wolketal 



20021) and LMC (IMaddox et all 120091) . IC131-se would simply be extreme in that the extent and 
amount of the non-thermal emission would be much larger than in any of the SNRs and Hn 
regions mentioned above. There are at least two objects of roughly compar able size to IC131-se 
for which n on-thermal compo nents have been reported. These are 30 Dor C (ISmith & Wangll2004l) 
and N51D (ICooper et al.ll2004l) w hose extents are ~ 90 x 100 pc and 120 x 200 pc, respectively. In 
30 Dor C, Smith & Wang " (|2004 ) found that roughly half of the total emission is thermal and half 
is non-thermal. The non-thermal emission is brightest along the western edge of the bubble and 
has a luminosity of about three times the total X-ray luminosity of IC 13 1-se. In N51D, about 30% 
of the total X-ray luminosity is non-thermal (ICooper et al.1 12004)) and has a luminosity of about 
one eighth of that observed in IC 13 1-se. The photon index that is requi red to explain the X -ray 
emission inIC131-se of T ~ 2.1 is consistent with the results reported by IBamba et al.l (12004) and 



Smith & Wang (2004) for 30 Dor C. 



Recent radio observations of M33 at 20 cm and 3 .6 cm revealed that about 7 - 30% of the radio 



emis sion in IC131 is synchrotron emission (0.29 mJy and 1.33 mJy, respectively. ITabatabaei et al 



2007|) . This radiation indicates the presence of relativistic CRs, which can also be the cause of 
synchrotron emission in the X-ray regime. Alternatively, non-thermal X-rays can also be produced 
by inverse Compton scattering of relativistic particles which interact with low energy photons from 
IC131 or the Cosmic Microwave Background. The challenge associated with a non-thermal inter- 
pretation of the emission we observe in IC13 1-se is that there are no detailed theoretical models for 
confronting the observations of IC 13 1-se or any of the other objects in which harder X-ray com- 
ponents have been observed. The basic idea is that multiple shocks from SNe and stellar winds 
reverberating through the bubble provide a fertile region for the acceleration of electrons and per- 
haps cosmic rays. It is generally believed that synchrotron radiation accounts for the non-thermal 
emission from SNRs. 

To determine which of the non-thermal mechanisms might be the dominant one, we start with 
an estimate of the radiative losses due to inverse Compton scattering. The total energy lost by one 
electron per unit time that passes through a photon field of energy density U rad can be written as: 



dE\ 
dt J ic 



-a T c(3 7 U rad 



(3) 



where <7t is the Thom pson cross-section, (3 = v/c, and 7 is the Lorentz factor (see, e.g jRybicki & Lightman 
19861 : LLongair Ill994|) . The stellar radiation field can be approximated by U ra( i = L stars /AncR 2 , 
where L stars is the luminosity of the stars and R is the radius of the X-ray emitting volume in cm. 
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Because an accurate measure of L stars is lacking, we assume a conservative upper limit for the 16 
B-type stars inside the bubble of 4.9 x 10 37 erg s -1 , which results in U raa - = 1.4 x 1CT 15 ergcm~ 3 . 
This energy density is 300 times smaller t han the energy density for the cosmic microwave back- 



ground of U C mb = 4.2 x 10 13 erg cm 3 (|Longair Ill994|) 



Because the relation for radiative losses from synchrotron emission is identical to Eq. ©, 
except that the radiation energy density U ra d is replaced by the energy density of the magnetic field 
Umag = B 2 /8n, the ratio of the two loss rates is simply (dE/dt)i C /(dE/dt) syn = U rad /U mag . 
Assuming a magnetic field strength across IC131-se of 20/xG (F. Tabatabaei, priv. comm.), the 
corresponding magnetic energy density is U mag = 1.6 x 1CT 11 erg cm -3 , yielding U mag ^> U ra d. 
Hence, synchrotron losses clearly dominate over losses from inverse Compton scattering for the 
adopted magnetic field strength and stellar radiation field. 

If we assume that the observed X-ray emission is synchrotron radiation, we can estimate the 
total power of the relativistic electron population {(3 ~ 1). We assume further that the number 
density distribution of the electrons follows a piecewise power law distribution (broken power 
law), where each piece has the form N("f) = iV 7~ p and extrapolate one part from the radio and 
the other part from the X-ray regime. Here, p is the power law index of the energy distribution of 
the radio and X-ray electrons, p is related to the radio photon index T r via p = 1 — 2T r , whereas 
p — 1 + 2T X holds for X-ray energies. For the radio, we adopt a photon index of T r = 0.3 
(Tabatabaei, priv. comm.) and for the X-rays we use T x = 2.13 (Table [T]). After replacing U ra d 
with U mag in Eq. ©, multiplying with N(*y), and integrating over 7, we obtain: 

P syn = 2.66 x 10- 14 U mag N ( 7 ™^~J mi " ) erg s'W 3 . (4) 

From Eq. ©, we can determine the normalization constants N for both parts of the piecewise 
power law distribution by integrating over the appropriate range of 7 for the radio and X- ray emis 



sion, assuming that the electrons emit at the critical frequency v c = 7 2 e£? /27i m P (e.g., Longair 
1994). The luminosities in these bands amount to L rad = 2.8 x 10 33 erg s _1 (iTabatabaei et al 



20071) and L x = 8.7 x 10 35 erg s -1 . 

The total energy, E tot , stored in the relativistic electron population can be written as: 



Etot = J N^E^djdV = N m e c 2 V , (5) 

where the first power law component is integrated from 7 min = to 'jbreak an d the second segment 
from 'jbreak to ^max = oo. The Lorentz factor at the break point of the broken power law can be esti- 
mated by extrapolating the power law segments to where they intersect, yielding ^break = 1.6 x 10 6 . 
The total energy in the synchrotron emitting electron population is therefore 2.83 x 10 51 erg which 
is ~14%-30% of the total energy provided by the estimated number of SNe (9-21). It follows 
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that non-thermal synchrotron emission can be a significant contributor to the overall emission. 
This result is somewhat unexpected and the challenge to this interpretation is to determine how 
a relatively strong magnetic field can be maintained on large sca l es in the absence of colliding 
winds from O stars. This result is also in contrast to iMuno et al.1 ([2006) who found that inverse 
Compton scattering is the most likely loss mechanism to explain the non-thermal X-ray emission 
in Westerlund 1 . 



3.3.3. Conclusions 



From an observational perspective, we have shown that IC131-se is an exceptional object in 
terms of the combination of the extent and the hardness of its X-ray spectrum. We have explored 
the feasibility of various emission mechanisms, and we conclude that, with the available data, we 
cannot distinguish between interpretations based on thermal emission from a hot ISM filling the 
interior of the superbubble or non-thermal radiation which seems to originate from synchrotron 
photons. The primary challenge set down by the presented data is for others to attempt to produce 
models of superbubbles and GHRs that yield high enough electron temperatures to reproduce the 
spectrum or enough high energy particles to generate the X-rays that are observed via synchrotron 
radiation. It would be interesting to know what makes IC131-se so peculiar, and if there are any 
other GHRs which are really similar. 



4. Point source CXO J013315.10+304453.0 
4.1. Imaging results 

The distribution of the diffuse X-ray emission in IC131-se in the soft, medium, and hard 
energy band (Fig.QJ) reveals several localized enhancements of the X-ray emission, some of which 
can be associated with optical counterparts. Based on the full ChASeM33 source catalog (which 
will be presented elsewhere), the only X-ray point source in IC131-se, whose extent is consistent 
with the PSF, is CXO J0133 15. 10+304453.0 (FL073). 

This object is visible in Fig. \Vp as a white circular region (see red plus sign) and seems to 
have an optical counterpart (Fig. [5]) which can be ide ntified as J0133 15. 17 +304452.5, a faint star 
of vis ual magnitude m v = 22.2 and B — V — 0.14 (hvlassev et al. Ibood). Using E(B - V) = 
0.07 ( Massey et al. 2006 ) and a distance of 817 kpc (|Freedman et al. I 200 1|) the apparent visual 



magnitude translates into an absolute magnitude of My = —2.61. Although optical spectroscopic 
data, which could clarify the spectral type and a possible membership to IC131-se is missing, 
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Fig. 5. — A zoom into the region highlighted in Fig. QJ revealing a possible optical counterpart to 
FL073. Contours are from the unbinned merged event list in the 0.35 -8.0 keV energy band and are 
plotted on a square root scale between 2 and 8 counts per pixel. See the electronic edition of the 
Journal for a color version of this figure. 



color and visual mag nitude seem to imply a star of spectral type B2V with a mass of about 9M Q 
(|Schmidt-Kalerj 1 1982b . Such high-m ass stars are often fo und as companion stars in high-mass 
X-ray binary (HMXRB) systems (e.g JMcBride et alJl2008h . 



4.2. Spectroscopic Results 

Because the point source has only ~150 counts, we used only absorbed single component 
models to fit the X-ray spectrum of FL073. The H I absorption was again modeled with a combi- 
nation of tbabs and tbvarabs models, assuming the same values for the Galactic and internal 
absorption in M33 as in Section |3~2l In Table[2]the best-fit results are shown for a non-thermal and 
a number of thermal plasma models. Fig. [6] presents the spectrum and 90% confidence contours 
of the adopted power law model. Because the point source appears to be embedded in diffuse 
gas, the background was taken from regions of the diffuse emission. As a consequence, a rela- 
tively hard background spectrum was subtracted which apparently softened the net emission and 
required only a marginal internal column density. However, the upper limit on the absorption is in 
good agreement with the total integrated column density. 

Unfortunately, the fit statistics and unabsorbed luminosities are very similar for all models and 
do not favor one model over another. Because there are no indications of emission line features and 
the spectrum looks relatively smooth, we assume the absorbed power law model to be the best fit 
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Table 2. Fit parameters (Point source) 



tbabs x (tbvarabs. . . ) a 


(pow) 


(brems) 


(apec) 


(nei) 


A, red 


1.09 


1.13 


1.11 


1.10 


degrees of freedom 


15 


15 


15 


14 


P value 


0.53 


0.58 


0.59 


0.52 


K (10~ 6 cm" 5 ) b 


1 20+ 048 


1 4Q+0.42 


4 46+ ' 79 

^■^°-0.74 


4 40+ - 73 

^• HU -0.75 


N H (10 22 cm" 2 ) 


<00 i+0.26 
^ u>ui -0 .01 


<0 01+ 017 


<o.oii^I 


<0.01_°; ? 


r c 


i OO+0.42 
L -°°-Q.26 








kT e (keV) d 




3 SQ+4.07 


3.961 3 ; 77 


3.981?-^ 


r (10 13 s/cm 3 ) c 








4 76 +4 - 76 

H - /D -4.76 


Fx (10- 15 erg/s/cm 2 ) f 


e o 7 +0.66 
3 - 8/ -2.03 


5 26 +0 - 71 


5 36 +0 - 57 


<r o 7 +0.92 
J - J 1 -2.52 


F x (10" 15 erg/s/cm 2 ) f 


6.411^1 


5.66_ 4>00 


5 74+0.6I 


5 74+0.98 
J - '^-2.69 


Lo.35-8 (10 35 erg/s)s 


i on+0.14 

l.Z8_ 044 


i 10+O.I8 
1 " 1J -0.8Q 




u4i°S2 


L 2 _ 8 (10 35 erg/s) h 


62+ 06 


55+ - 07 

u - JJ -0.30 


o.58i°:« 7 


0.5818S 



a See Table Q] for a description of the parameters. 



IC131 point source tbabs(tbvarabs(pow)) Confidence Contours (90%): tbabs(tbvarabs(pow)) 
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Energy (keV) Parameter: nH (1 22 cm- 2 ) 



Fig. 6.— Shown are the spectrum of CXO J013315. 10+304453.0 (FL073) and the best model fit 
of the absorbed power law (left panel) together with the 90% confidence contours (right panel). 
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Fig. 7. — Long-term variability of FL073 (the inlay shows the full lightcurve, including the archival 
ObsID 1730). 



to the data. The unabsorbed broadband X-ray luminosity of FL073 amounts to about 1.3 x 10 35 erg 
s~\ ~50% of which are emitted at hard X-rays above 2.0 keV. 



4.3. Discussion 



The photon index of FL073 of T — 1. 881°' te, is consi s tent with that of two transient X-ray 



sources in M33, called XRT-2 or XRT-5 in IWilliams et al.l (I2008|) which have photon indices of 
1.6±0.4 and 1.9±0.5, respectively. XRT-2 seems to have an optical counterpart which could be an 
early B star with very similar color (B — V = 0.1) and visual magnitude (my = 22.8) compared 
to the counterpart of FL073. The XRB hypothesis is corrobo rated by FL073's hardness ratios of 
HR1 = —0.17 and HR2 = —0.13 (see IPlucinsky et al.1120081 for definition and adopted energy 
bands). They agree well with those of the other XRB candidates that can be fitted with an absorbed 
power law < 2.0 and an Nu of ~ 10 21 cm -2 (see their Fig. 5). Although the photon index and the 
hardness ratios for FL073 could also indicate an AGN, the column density suggests that this object 
is part of IC131-se and is likely not a background AGN. 



veys of M33 (IPietsch et all 12004 



Based on the data presented b y the ROSAT (IHaberl & Pietschl 1200 lh and XMM-Newton sur- 



Misanovic et al.l 



2006h it cannot conclusively be determined 



whether FL073 was detected, as the poor spatial resolution of these satellites did not allow for 
a clear separation of poin t sources and diffuse emission. FL073 is not listed in the M33 CFHT 
optical variability survey (IHartman et al.1 120061) as a variable source. In X-rays this source is de- 
tectable in all ObsIDs and reveals a possible long-term variability (Fig. [7]). Fitting the lightcurve 
with an assumed constant flux yields xV* 7 = 9.47/3, which corresponds to a 3% probability of 
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being constant. Although the amplitude of the variation is not strong (~ 2.2a), we suggest that 
FL073 is associated with a HMXRB system in which the optical companion is a B2-type star with 
a mass of ~9M Q . It should be noted that only the high angular resolution of Chandra allowed a 
proper treatment of this source whose luminosity is ~16% of the total diffuse X-ray luminosity. 
At lower resolution this point source would have seriously contaminated the emission from IC131. 



5. Summary and conclusions 

Compared to known superbubbles and Hn regions IC131-se appears to be peculiar with re- 
spect to its lopsided X-ray morphology, the large linear extent of the X-ray emission, the lack of 
massive O stars, its high electron temperature (in case the X-ray gas is thermal), and large fraction 
of hard X-rays. The X-ray spectrum of the extended emission in IC131-se can be equally well 
fitted by an absorbed power law (T ~ 2.1) or an absorbed thermal plasma model (kT e ~ 4.3 keV). 
These models predict a total (0.35- 8.0 keV) unabsorbed X-ray luminosity of about 9 x 10 35 erg 
s _1 and 7 x 10 35 erg s~\ respectively, with 39% and 53% of the luminosity being emitted above 
2keV. Apparently IC131 possesses not only the hardest X-ray spectrum among the known GHRs, 
but, in case the gas is thermal, also the hottest X-ray plasma. 

Besides the facts that thermal X-ray emission is expected from GHRs and that there might 
be weak emission lines below 1 keV, the disagreement between the number of SNe progenitors 
obtained from the IMF (9-21) and that derived from the internal thermal energy ( 1 1 8±59) argues 
against a pure thermal origin of the gas. We estimate a low electron density of n e ~ 0.05 cm -3 
(assuming f x = 0.9) and an X-ray gas mass of about 3300 M . The age of the stellar population 
inside the bubble seems to range from 8 to 77 Myr and given a cooling timescale of about 1 Gyr, 
the gas did not have time to cool significantly. If the gas is thermal, the standard bubble models 
require shocks from SNRs and SNe to heat the gas. The fundamental problem which argues against 
a thermal plasma is that we are not aware of a mechanism for GHRs which can produce electron 
temperatures anywhere close to kT e = 4 keV. 

Non-thermal X-ray emission is also known to originate from SNRs, GHRs or superbubbles 
and the spectral power law fit is statistically as good as the one for a thermal plasma. Synchrotron 
emission as opposed to inverse Compton scattering is a plausible mechanism for the assumed mag- 
netic field strength of B = 20fiG and the total energy of the relativistic synchrotron-emitting elec- 
tron population is about 14% to 30% of the total energy provided by SNe. If a non-thermal com- 
ponent exists, synchrotron losses clearly dominate over losses from inverse Compton scattering. 
For a purely non-thermal origin, however, it remains to be investigated what kinds of mechanisms 
are able to produce non-thermal X-rays or accelerate particles to relativistic energies on scales of 
200 pc. To allow for a more quantitative analysis of the non-thermal X-ray emission, models are 
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required which predict the fraction of non-thermal radiation from GHRs. A combined thermal and 
non-thermal model seems to be generally appropriate, too, irrespective of the unreasonably high 
column density and X-ray luminosity derived from those fits. 

Unfortunately, with the present data the nature of the extended X-ray emission cannot be con- 
clusively determined. Clearly, IC131-se is an important object and challenges the standard bubble 
model as well as our understanding of CR acceleration in superbubbles. It is remarkable and re- 
mains to be understood why objects similar to IC131-se have not been detected before. Future 
investigations would greatly profit from more sensitive and deeper X-ray observations (as for ex- 
ample provided by IXO), high-resolution optical spectrophotometric data, and models which can 
explain electron temperatures > 1 keV and make predictions on the non-thermal X-ray emission 
in GHRs. The deeper X-ray observations could establish the presence of X-ray emission lines and 
settle the question on the nature of the emission, while a detailed stellar population analysis could 
provide a more accurate IMF. Clearly, NGC604 and IC131-se seem to be in completely different 
evolutionary stages. All O-type stars in IC131-se seem to have exploded as SNe, whereas the 
western part of NGC604 awaits the first SNe to occur. 

We detect only one X-ray point source in IC131. This source, CXO J013315. 10+304453.0 
(FL073), appears to be time variable and is possibly a HMXRB with an optical counterpart which 
could be a B2V star with a mass of about 9M Q . 



This work has made use of SAOImage DS9 (|Joye & Mandel 1120031) . developed by the SAO, 
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